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1.  SUMMARY 


The  principal  objective s  of  this  work  cope: 

(1)  to  inveetigate  ray  theory ,  normal  mode  or  my  newly-developed 
method  by  comparing  the  computed  reeulte  with  field  measured 
data  for  predicting  underwater  sound  propagation  lose  in  the 
Arctic ,  and 

(2)  to  apply  appropriate  taahniqusz  to  evaluate  acoustic  intensity 

and  transmission  lose  at  promising  locations  from  a  nmber  of 
sound  souross  in  the  Arctic.  .  ,  ,  *  » 

a  r, .  y 

Because  of  the  limited  time  available  for  completion  of  the  project, 
efforts  to  fulfil  both  objectives  had  to  be  pursued  at  the  same  time. 
Thus  while  a  program  of  research  using  direct  integration  (FFP,  Sec¬ 
tion  6)  was  being  carried  out  at  Lamont-Doherty  Geological  Observatory, 
concurrently  a  ray  tracing  program  geared  specifically  for  the  calcula¬ 
tion  of  acoustic  intensities  in  the  Arctic  was  being  implemented  at 
Aerophysiae  Research  Corporation. 

Details  of  the  direct  integration  method  are  given  in  a  separate  report 
(Reference  5,  Section  6).  Comparisons  of  results  calculated  by  ray 
theory  and  by  direct  integration  are  illustrated  on  page  89. 

The  Arctic  Ocean  bottom  is  graphically  illustrated  on  pages  61-65  by 
plotting  the  various  depth  contours.  These  and  other  considerations 
have  led  to  the  choice  of  a  nmber  of  sound  souross  and  receivers, 
which  are  partitioned  into  five  groups  roughly  in  accordance  with  the 
bottom  topography.  The  location  of  the  sources  and  receivers  is 
tabulated  and  indicated  in  the  contour  maps  on  pages  87-78. 

Comparisons  of  field  measured  loss-range  curves  with  calculated  results 
by  ray  theory  are  shown  an  pages  77-88.  The  complete  computer  output 
of  acoustic  intensities  and  transmission  loss  is  presented  in  Part  II 
of  the  report.  It  is  indeed  gratifying  to  note  that  on  the  thole 
the  computed  results  agree  well  with  field  data. 


\0J 

I  ** 


As  follow  up  work  in  the  future  the  following  suggestions  will  he  made. 

(1)  Several  refinements  can  be  experimented  and  explored  in  the  ray 
tracing  program.  Correction  for  acoustic  intensity  for  regions 
close  to  the  caustic  zones  should  he  re-examined.  Also  the  bot¬ 
tom  reflection  loss  calculations  can  be  reformulated  by  utilizing 
more  information  of  the  bottom  topography.  Similarly  the  direct 
integration  method  can  be  developed  further  to  take  into  consi¬ 
deration  a  varying  bottom. 

(2)  Site  survey  (propagation  loss  calculations)  will  be  carried  out 
at  various  frequencies  in  the  band  10  to  500  hz.  Specifically 
computations  are  to  be  performed  by  both  ray  theory  and  direct 
integration  as  cross-check  at  10,  15,  20,  30,  50,  60,  80,  100, 
150,  200,  250,  300,  400  and  500  hz. 

(3)  From  the  results  of  (2)  the  detection  ranges  for  typical  surface 
and  underwater  sources  will  be  estimated  from  the  sonar  equation. 
Detection  systems  which  appear  feasible  for  use  in  the  Arctic 
will  be  considered,  viz.,  various  types  of  hydrophone  arrays, 
signal  processing,  data  display  and  data  gathering. 

The  fact  that  now  rapid  and  efficient  prediction  can  be  made  on  the 
underwater  sound  transmission  loss  indeed  paves  the  way  for  future 
exploration  and  research  in  the  Arctic  Ocean  area. 


2.  THEORETICAL  DEVELOPMENTS 


2.1  ENERGY  ALONG  A  RAY  PATH 

Consider  two-dimensional  ray  propagation  in  the  (x,y)  plane  with  sound 
velocity  c  =  c(x,y) 


y 


Fig.  1 

Let  s  be  the  distance  measured  along  a  ray  and  let  n  denote  the  local 
orthogonal  direction  as  shown.  If  0  is  the  slope  angle  of  the  ray,  the 
(.s,  nj  vectors  are  obtained  by  rotating  the  (x,  y)  vectors  counterclock¬ 
wise  through  0. 


The  curvature  1/R  is  given  by 


I 

R  ds 

Since  the  ray  equations  read 

df*  „  fc  dx  dfc.  dx 

^c2  ^  x  ds  y  ds^  ds 

(c 


ds2 

d2y 


x  ds 
dx 


ds2 


'x  ds  T  Kv  — ) 


-  C, 


-  C, 


CD 


(2) 
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It  follows  easily  that 


c 


de 

ds 


cx  sin0  -  Cy.  cos0  = 


so  that 


1  3  , 

R  *  to  Wn  c)  C3) 

The  rays  are  orthogonal  to  surfaces  of  constant  phase;  for  brevity,  such 
surfaces  will  be  termed  wave  fronts.  Denote  the  local  curvature  of  a 
wave  front  by  (1/RW);  taken  as  positive  n  the  rays  are  diverging.  The 
immediate  purpose  is  to  compute  d/ds(l/Rw). 

For  convenience,  introduce  an  orthogonal  curvilinear  coordinate  system 
with  £  constant  on  a  ray  and  n  constant  on  a  wave  front. 


Clearly, 


Thus, 


cos0  r  .  „  Sr-yrnCose 

xny5  ?ncose  '  —  -  056r>sinej  (6) 


Similarly, 


d  -  C°S0  To  «  0nXrnCOS0 

dn  V  -  xny5  t05ncos0  £ -  -  e^sino]- 


combining,  we  obtain 


)  =  d_  f  K  cos0  ,xgn  ynu 

ds  Rw  dn  x  y  ^  R  “  R  ^ 


Next,  differentiate  the  first  of  Equations  (4)  with  respect  to  £  and 
the  second  with  respect  to  n  ;  then  solve  the  resulting  equations  for 
y£n  and  x£n  t0  obtain  (with  the  help  of  Equations  (5)) 

y»  =  rL  .  tanO> 

COS0  v,rw  r  > 


x  =  Xr,yE  r  i_  tanO. 
5n  cose  l"r  "  rw  J 


Substitution  into  Equation  (8)  now  gives  the  desired  result: 


1_  (L.)  =  1_  .1  _  1_  i 

ds  VRW7  dnV  R2  '  R2 


To  obtain  a  convenient  formula  for  d/dn(l/R),  we  observe  that  Equation  (3) 
yields  the  following. 


c)x  sine  -  {In  c)  cose] 


=  [(Jin  c)  cose  +  {In  c)  sine] 


de 

dn 


"  [(^n  c)  sin2e  -  2  (Jin  c)  sine  cose  +  (Jin  c)  cos2e] 
aa  xy  w  ■ 


yy 


=  ^  [zran  c)]  ■  ci 


(ii) 


where  d2/dn2  is  defined  in  the  "straight  line"  sense  via  the  progression 
between  the  last  two  lines  of  Equation  11.  Carrying  out  the  differentiation 
of  (Jin  c)  gives  the  alternative  form 


—  (— ) 
ds  V 


1  dc 
cRw  ds 


1_  d2c  1 
c  dn2  "  Ry2 


(12) 


where 


dH?  =  cxx  sin2fl  "  2cxy  sine  cose  +  cyy  cos2e  (13) 

Suppose  now  that  the  source  is  located  at  some  point  on  the  (negative) 
y  axis.  At  each  point  along  a  ray  we  know  the  range  x,  the  slope  angle 
0,  and  the  wave  front  curvature  1/RW  (as  a  result  of  Equation  (12)). 

Let  F  denote  the  intensity  (energy  rate  per  unit  area)  along  a  ray.  It 
now  follows  from  simple  geometry  that 


.1 

F  ds 


cos6  1 
x  Rw 


(14) 


If  A(x,  y)  denotes  the  acoustic  attenuation  in  the  water,  then  the 
final  equation  becomes 


1  =  £££!  +  1_  +  A(x  v) 

F  ds  x  Rw  1  ,yJ 


(15) 


In  a  program,  it  is  worthwhile  to  record  each  term  in  Equation  (15) 
separately;  thus,  F  is  divided  into  the  three  terms  Fr,  F^,  and  F^ 
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via 


F  •  Fr  •  FW  •  Fa 

so  that 

1  dFR  _  cos6  ^  1  dx 
Fr  ds  •  x  x  ds 


1  '  dFvv  _  1_ 

'FW  ds  “  Rw 

1  dFA  « r  , 

—k  ar  ■  A(x-  ^ 

The  first  of  these  equations  yields 
-  const . 

r  =  — 

R  x 


(16) 


when  the  constant  is  chosen  to  be  the  intensity  of  the  source  at  unit 
distance  (i.e.,  s  =  1),  along  the  chosen  ray;  thus,  Fw  and  FA  are  each 
unity  at  that  point. 


A  knowledge  of  F^  has  an  interesting  physical  interpretation.  Consider 
two  rays  emanating  from  the  source  at  angle  0q,  at  an  incremental  angle 
d0Q  apart.  At  unit  distance  from  the  source,  F^  *  1.  At  a  terminal 
point  xt,  where  the  slope  angle  is  0t  and  the  value  of  Fw  is  F^t*  we  can 
use  the  fact  that  F^  is  inversely  proportional  to  the  normal  spacing  dn 
between  adjacent  rays  to  write 

(l)(l)d0o  =  FWt|dn|  (17) 


so  that  | dn/d0Q |  =  1/F^t 
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2.2  CORRECTIONS  TO  RAY  ACOUSTICS 


Let  (xj,  yi ,  zj)  be  rectangular  Cartesian  coordinates  and  t  be  time.  Let 
SL  be  some  characteristic  length  (e.g. ,  water  depth),  and  define  the 
non-dimensional  coordinates  by  x  *  Xj/L,  etc.  Consider  the  wave  equation 

•  ■  4  *tt  08) 
cz 

where  <j>  is  the  velocity  potential,  c  -  c(x,  y,  z) ,  a  function  of  the  space 
coordinates.  Let  cQ  be  some  constant  which  can  be  taken  approximately  as 
the  average  value  of  c  over  the  region  of  interest. 

For  periodic  motion  we  assume 

$  =  Re  (u  e*a,t)  (19) 

whore  uj  is  the  (constant)  angular  frequency;  u(x,  y,  z)  is  the  complex 
amplitude.  Substitution  of  Equation  (19)  into  (18)  shows  that  u  satisfies 
the  equation 


Au  +  k2n2u  ■  0 


(20) 


where 


lL 

8y2 


and 


n 


c 


The  quantity  k  is  non-dimensional  and  in  practice  k2  is  usually  »1.  For 
example,  at  a  frequency  of  20  Hz  and  with  a  depth  of  2,000  feet,  k2  is 
approximately  2500. 


In  acoustic  propagation  problems  one  is  often  dealing  with  a  sequence  of 
progressing  waves  of  approximately  constant  phase;  for  example,  plane  waves, 
u  *  ei^x,  and  spherical  waves,  u  ■  e“i^r/r.  For  such  problems  a  natural 
substitution  is 


u  ■  a(x,  y,  z)e 


ik8(x,y,z) 


(21) 


I 

I 

where  the  function  3(x,  y,  z)  is  to  provide  the  travelling  wave  character. 
Substitution  of  Equation  (21)  into  (20)  gives 

k2[n2-(Sx2+Sy2+Sz2)]a  +  ik[2(«x3x+o.y{3y+ozBz)  +  ahfi]  +  Aa  =  0  (22) 

For  large  k  the  first  term  is  dominant  so  that  as  a  first  approximation 

ex2  *  Sy2  +  Sz2  =  "2  (23) 

This  is  the  eikonal  equation,  and  its  characteristics  are  the  acoustic  rays. 
If  s  is  the  distance  along  a  characteristic,  the  differential  equations  of  a 
characteristic  become 


4*  =  ix  ,  3*  dz  _  3j 

ds  n  ds  n  1  ds  n 


dSx 

3s~  =  n*' 


dgz 

ny  z r =  nz 


from  which  the  conventional  ray  equation  follows: 


d2x  1  dx  ,  dx  A  dy  dz^  1 

ds2  a  ds  (nxds  Vds  +  ^d?5  +  n  n> 


(24) 


(25) 


The  term  that  is  next  important  in  (22)  can  also  be  put  to  zero  to  get  a 
second  approximation  by  choosing  o  such  that 

2(axSx  +  ayS y  +  oz8z)  ♦  aAS  «  0  (26) 

From  (24)  and  (26)  it  follows  along  a  ray 

2n~  +  oAB  =  0  (27) 

div  (a2  gTad  S)  =  0  (28) 
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and  by  the  divergence  theorem  for  a  closed  surface  S 

I  a2  grad  3  *  dS  =  0  (29) 

S 

Consider  a  section  of  a  ray  tube  and  make  use  of  (29) .  The  contribution  to 
the  integral  from  the  tube  sides  is  zero  since  grad  3  is  in  the  direction 
of  the  ray;  so  that  if  the  two  end  sections  of  the  ray  tube  be  given  the 
indices  1  and  2,  we  have 

I  a2ndS  =  /  a2ndS  •  (30) 

Si  S2 

Since  the  rate  of  the  energy  propagation  per  unit  area  is  proportional  to 
(amplitude)  /c,  (30)  implies  that  the  rate  of  energy  propagation  is  constant 
along  a  ray  tube. 


This  leads  to  the  conclusion  that  the  use  of  ray  theory  is  equivalent  to  the 
omission  of  the  last  term  in  (22).  One  way  of  examining  the  adequacy  of  the 
ray  theory  is  to  solve  the  problem  by  ray  theory  to  obtain  a  and  3  at  each 
point  and  to  examine  the  magnitude  of  Aa.  Another  approach  is  to  iterate, 
compute  Aaa  f,  say,  and  recompute  a  by 

2n^|  +  aA3  +  ~  f  =  0  (31) 

A  comparison  of  the  revised  values  of  a  should  provide  a  criterion  for  the 
adequacy  of  the  ray  approximation. 


Such  iterations  can  be  carried  out  more  systematically  by  using  the  asymp¬ 
totic  expansion  below. 
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Substitution  into  (22)  with  S  still  satisfying  (23)  yields  the  sequence 


of  equations: 

,  daC°> 
2nds 

+  ct^AS  =  0 

2n  ““(1) 
2n  ds 

+  a^AS  +  Aa^  =  0 

2n  f t2) 

+  ct^AS  +  Aa^1^  =  0 

(33) 

so  that  , 

«(1>,  .  . 

may  be  determined  in  succession.- 

Note  that  the 

ray  paths  are  unaffected  by  this  iterative  procedure. 
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2.3  NORMAL  MODES 


i 


p 

» 


■ 


> 


UTe  summarize  here  some  of  the  relations  between  the  methods  of  ray  tracing, 
normal  modes ,  and  image  methods . 

Consider  first  the  problem  of  an  infinite  ocean  of  constant  depth,  in  which 
c  -  c(z).  Let  there  be  a  concentrated  source,  of  angular  frequency  oi,  at 
a  distance  d  below  the  surface,  and  introduce  a  coordinate  system  as  shown: 


1 


surface 


bottom 


Fig.  3 

Writing  $  -  i|i(r,z)e1<ut,  the  governing  equation  becomes 
*rr  *  ~  *r  *  *zz  *  ^  *  ■  lfe:eL^Z'dJ 


(34) 


where  e  >  0,e  -*■  0.  Here  6(x)  is  the  delta  function.  As  boundary  conditions 
we  will  take  the  special  situation: 


Kr,  0)  =  0, 

>i'z(r,  H)  =  0,  (35) 

in  order  to  simplify  comparisons  between  the  various  methods. 
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(1)  Normal  Mode  Method 


We  write  <Hr»  z)  =  £an(r)  fn(z)  where  fn(z)  are  appropriate  expansion 
functions  chosen  by  the  use  of  the  associate  Sturm- Li ouvi lie  problem: 


I  f  +  C^2~  '  ~~  +  X)  f  =  0 


i 


f(Q'j  =  0, 


f'(H)  =  0. 


(36) 


Here  eg  is  a  constant  equal  to  some  average  value  of  c(z) .  Denote  the 
eigenvalues  by  Xn,  n  =  1,  2,  .  .  . ,  and  the  corresponding  eigenfunctions 
by  fn(z)-  Note  that  if  c(z)  ■  eg,  we  have  Xn  =  (2n-l)2u2/(2H)2  and 
fn  =  sin(X„  z) .  The  orthogonality  condition  requires 
H 

/  fm  fn  dz  «  0  for  Xm  /  Xn 


(37) 


Let 


H 


/  fn2  dZ  a  gn. 


(38) 


then 


H 


an(r)  -  ~  /  1»(r,  z)fn(z)dz 


En 


(39) 


Multiplying  Equation  (34)  by  fn(z)/gn  and  integrating,  we  obtain  by  integration 
by  parts 


*'"*?*;*' ft 


0  "  "  2irc  «n 

If  Xn  <  u2/c20  (say  for  n  <  N) ,  the  solution  is  a  multiple  of 

Jo(r  ™)2/c20  -  Xn)  for  r  <  e 


(40) 
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and  of  (r/u2/c2Q  -  A^  for  r  >  e;.  u 

t  I 

‘  •  '  *  *  »  *  »•  W  I  t  I'  »  •  I  till. 

matching  the  condition  at  r  ■  e,  and  permitting  e  -*-0  leads  to 

an(r)  “  "  4?&r  H^2)  (r*^2/c20  "  V  for  n  <  N  (41) 

Here  N  is  the  upper  limit  of  values  of  n  for  which  u2/c20  >  xn.  Similarly, 
for  n  >  N,  we  obtain 

an(r)  "  "  K0(r^Xn  -  “2/c20  n  >  N  (42) 

For  large  r,  Kq  is  exponentially  small  so  that  these  modes  are  relatively 

unimportant. 

As  a  numerical  example,  consider  again  c  »  Cq  ■  5000  ft/sec.,  say,  and  take 
u/c0  =  2ir(lo)/5000  -  2ir/500.  Then  An  -  [ (2n-l) w/(2H)  J 2  .  Let  H  »  1000  feet 
then  fn  *  sini^  z,  gn  «  H/2,  and 

<J»  (r.z)  ■  I  C-  sin/xjj  d)  (r/«o2/c2Q  -  A„)  sin/A,,  z 

+  ^  wH  sin^An  sin^\i  (43) 

5 

The  ratio  of  mode  5  to  4  is  of  the  order 

lcot25oo  ^  /  ho(2)  '5500  **  r>  *  C44) 

for  large  r,  so  that  for  r  beyond  1000  feet,  say,  we  can  neglect  all  modes 
beyond  the  fourth.  For  a  frequency  of  1  cps  in  this  example,  all  modes 
decay  exponentially;  for  a  frequency  of  1000  cps,  400  modes  are  required. 


* 


We  choose  H^2^  rather  than  in  order  to  obtain  asymptotically  an 
outgoing  wave  for  large  r,  in  conjunction  with  the  time  factor  eiwt'. 


(2)  Images  or  Rays 


A  delta  function  source  at  (x0,  y0,  z0)  in  an  infinite  medium  satisfies  the 
equation  (for  c(z)  =  Cq  ) : 

I 

At  +  —  t  =  6(x  -  xQ)  6  (y  -  y0)  6  (z  -  zQ)  (45) 

C  0 

Let  p  denote  the  distance  from  this  source ;  th* ' 

imp 

*  ■  -  4^  0  C°  (46) 

(Even  for  c  =  c(z)  this  is  adequate  for  small  p). 

The  problem  considered  on  Page  11  can  be  interpreted  in  terms  of  coupled 
pairs  of  images;  for  the  case  c(z)  ■  Cq,  tho  analysis  is  particularly  easy: 


Ut 


It  follows  that 


4> 


1  r1 
7—  { —  e 
4tt  px 


e  c0  -1—  e  c° 
P3 


iwS2 


(47) 


This  solution  is  identical  with  that  obtained  by  the  ray  method  since  there 
is  perfect  reflection  at  top  and  bottom  with  a  phase  change  of^  at  the  top 
surface.  Besides,  the  following  relations  hold: 


Px2  ■  r2  +  (z  -  d)2 
P22  ■  r2  ♦  (H  -  z  ♦  H  -  d)2 

s?2  ■  r2  ♦  (z  ♦  d)2  (48) 

Advantage  can  be  taken  of  the  phase  differentials  to  simplify  the  summation 
of  Equation  (47) 


(3)  Hankel  Transform 

In  some  layer  problems,  preliminary  Fourier  decomposition  in  z-direction 
may  be  awkward.  An  alternative  method  is  that  of  Henkel  transform  which 

will  be  made  use  of  to  illustrate  this  problem.  Returning  to  Equation  (34) 
define 

00 

<Kp.z)  *  /r  J0(r^)\jj  (r,z)dr  (49) 

0 

and  (34)  then  implies 

*22  *  <7  -  e2)*  ■  h  s<2-’1)  '  (50) 

c 

which  is  an  ordinary  differential  equation  in  z. 

For  the  case  c(z)  «  cQ,  we  obtain 

p,z 3  ’  “  4ry  cosyH  {sinY(z*d-H)  ♦  siny(H  - 1 z-d | ) }  (51) 
where  y  «  /u2/cq2-p2  . 

In  any  event,  whether  c  -  c(z)  or  c  -  cq,  we  must  eventually  compute 

<»(r»  z)  ■  /  p  Jg(rp)  }(p,  z)dp  (52) 

0 

where  *,  in  (51)  may  well  have  singularities  on  the  path  of  integration 
(where  cosyH  ■  0). 

There  are  two  useful  devices  associated  with  (52).  The  first  of  those  involves 
the  addition  of  a  small  amount  of  damping  to  the  wave  equation  predecessor 
of  (34),  so  as  to  eventually  move  the  singularities  off  the  path  of  inte¬ 
gration.  This  results  in  the  term  (u2/c2)  in  (34)  becoming  replaced  by 
(<d2/c2  -  whore  c  is  a  small  positive  parameter;  eventually  e  ■*>  0.  The 

offect  of  this  in  the  special  case  (51)  is  to  locate  the  singularities  at 
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p  ■  -  *^2/c2  -  ei  -  ((2n-l)ir/2H)’2 


(53) 


for  integral  values  of  n.  These  are  a  finite  number  lying  just  below  the 
positive  real  axis,  an  infinite  number  lying  just  to  the  right  of  the  negative 
imaginary  axis  and  also  at  the  negatives  of  all  these  positions. 

The  second  device  is  to  recognize  that,  in  (52) ,  we  can  write 

J0(rp)  -  j  (H^frp)  +  Hq2^  (rp) },  (54) 

where  we  will  choose  the  branch  cut  of  the  Hankel  functions  to  lie  in  the 
upper  half  p-plane.  Then  as  we  move  continuously  from  the  positive  to  the 
negative  real  axis  through  the  negative  half  plane,  we  find  that  for  a  >  0, 

H(J2)(-ra)  -  -H0Cl)(ra), 
so  that  (52)  becomes 

*(r.z)  "  1  J  p  H0J)  (rp>  *•  (P»z)dp  ♦  }  J  p  H<2)  (rp)  «Kp,z)dp 
*  J  PH02)  ('rp)  *(p»z)dp  *jJ  PH0(2)  (rp)  KP,z)dp 

*  2  J  ^  ^rP)  KP»z)dp 

(2) 

where  Hq  has  no  branch  cut  in  the  lower  half  plane.  For  the  special  case 

of  (51)  a  straightforward  residue  calculation  (using  the  device  of  (S3)  and 

noting  that  the  point  p  -  u2/c2  is  not  a  branch  point  as  c  +  0)  reproduces 
the  results  of  method  (1). 


(4)  Perturbations 


An  efficient  way  in  which  to  handle  the  variable  c  case  via  normal  modes  would 
appear  to  be  by  use  of  perturbations.  Consider  the  mode  in  (36)  and 
write 


£  +  [et£h(zl  +  x  ]f  =  0> 

J  c  n  J  J 


where  c(z)  =  CQ//l+eh(z)  .  With  0  <  e  «1  we  also  write 

f, ♦  a®  ...  .. 

J »,  =  X<°>  *  cxj1)  ♦  e*X<2>  *  . 

so  as  to  obtain  the  sequence  of  equations 


£C°)„  +  XC0)  £(0)  = 

3  3  3 

far  +  xco)  £cd  a  .hfco)  .  X(D  £co) 

3  3  3  3  3  3 

f(2)"  +  x(0)  f(2)  _hf(l)  .  x(l)  f(l) 

j  J  3  j  3 


XC2)  f(0) 

3  3 


The  boundary  conditions  are 


fjp)  (0)  =  fjp)'  (H)  =0  for  all  p. 


(56) 


(57) 


(58) 


(59) 


It  follows  that  xj°>  =  (2j-l)2m2/4H2  and  that  fj0)  =  sin((2j-l)mz/2H) . 
Because  of  the  structure  of  the  second  of  (58) ,  there  can  be  a  solution 
only  if  its  right-handed  side  is  orthogonal  to  fj^\  thus  xj1)  must  be 
given  by 

xj13  *  -  |  /H  h(z)  sin2#^  z  dz  (60) 


K®  CUn  n°W  find  fj13*  discarding  any  multiple  of  fj0^  which  might  appear 
For  example,  let  h(z)  =  z,  y  =  (2j-l)ir/2H,  then 


and 


4r(l)  z2+2A 

'  4y 


Cl). 


cosyz  "  4^2  sinyz 


(61) 


(62) 


In  Equation  (41)  we  now  obtain  a  value  for  an(r)  ,  correct  to  orders,  by 
replacing  xn  by  ,  eXU),  %(d)  by  fC0)(d)  .  £fa)(d)  ^  by  revising 

fin  appropriately.  Higher  order  corrections  may  be  obtained'  similarly ; 
ot  course,  in  this  special  case  an  exact  analytical  solution  is  also  avail¬ 
able  but  our  main  purpose  was  to  outline  the  general  perturbation  technique 
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2.4  FOCAL  ZONES  AND  CAUSTICS 


In  the  neighborhood  of  a  focal  zone  or  caustic  some  distance  from  a  source 

in  a  cylindrically  symmetric  problem  the  motion  may  be  considered  two-dimensional. 


Consider  first  a  focusing  situation  in  which  a  ray  bundle  is  approaching  the 
origin.  We  write  4  ~  Re  (uelu)t)  with 

Au  +  k2u  =  0  (63) 

where  k  =  to/c;  we  take  c  as  constant  in  the  neighborhood  of  the  focus.  A 
plane  wave  solution  of  (63)  is  given  in  polar  coordinates  by 

u(r,  8)  =  elkr  “S<6-Y>  (64) 

where  y  is  the  angular  direction  from  which  the  wave  is  coming.  A  superpo¬ 
sition  of  such  waves  is  also  a  solution;  consider  therefore  a  wave  bundle  of 
the  form: 


u(r,  b)  =  /f eik ["*CY)-rcos ( e-y) ]  dy  (65) 

0 

where  f(y)  and  m(y)  are  arbitrary  real  functions.  Here  r  is  to  be  the  angle 
within  which  rays  approach  the  origin. 


A  useful  tool  in  dealing  with  integrals  like  that  in  (65)  is  the  method  of 
stationary  phase  which  states  that  if  there  is  one  point,  tg»  in  the  interval 
(a,  b)  at  which  h'(t)  =  0,  then  for  large  A  >  0, 


jh  eiAh(t) 


g(t)  dt  y*|h"(t0)  | 


g(tQ)e 


iXh(to)iiir/4 


(66) 


where  the 


+ 


sign  is  chosen  for  h"(tQ)  >  0  or  <  0,  respectively  (Reference  1). 


For  large  values  of  r,  Equation  (64)  thus  implies 

u(r»  6)  -J»  f(e)eiktm(6)+r]-i7r/4  (6?) 

for  0  <  6  <  r,  for  0  outsid>  this  range  (and  outside  the  range  0  <  0-ir  <  r) , 
we  have  u(r,  0)  -  0  (relatively).  Thus,  the  exact  solution  (65)  corresponds  to 
a  ray  bundle  of  aperture  r.  Let  us  choose  the  bundle  aperture  sufficiently 
small  that  f(e)  (the  amplitude  factor)  is  constant;  we  will  also  make  the 
slight  additional  simplification  m(0)  =  constant.  Then 

u(r,  0)  .  f  eikm  /eikr  Cos(0-Y)  dY  •  '  (68) 

0  1  1 

where  f  and  k  are  found  from  the  amplitude  and  phase  of  u  at  some  large  distance 
r,  as  given  by  (67). 


We  observe  first  that  on  the  other  side  or  the  focal  point,  the  method  of 
stationary  phase  can  also  be  applied  to  give,  for  large  r, 

u(r.  6)  -M  f  °ik(m'4M,r/4  (69 

The  change  in  sign  of  the  r  term  in  the  exponent  simply  corresponds  to  the 
fact  that  the  waves  are  now  outgoing  rather  than  incoming;  the  change  from 
-tt/4  to  7t/4  means  that  the  phase  has  advanced  by  ir/2  in  the  passage  through 
the  focus.  Observe,  however,  that  apart  from  these  effects  there  is  no 
amplitude  change;  the  energy  flow  rate  between  a  pair  of  converging  rays  is 
the  same  as  that  between  the  corresponding  pair  of  diverging  rays. 


Equation  (67)  is,  of  course,  not  valid  at  the  focal  point  itself— it  would 
give  u  =  ».  However,  (68)  is  valid  everywhere  so  that  at  r  »  0 


u 


f  eikm  r 


(70) 


Thus  the  ratio  u(0)/u(r),  r  large,  is  given  by 


*  /L_ 

a  u(r)  'i  2n 


(e-ikr»l./4)r 


(70) 


for  a  ray  bundle  of  aperture  r,  (the  factor  /r  simply  reflects  the  fact  that 
u(r)  -  l//r  as  we  move  along  a  ray).  Values  of  u  at  other  points  in  the 
immediate  neighborhood  of  the  focus  can  be  found  by  numerical  or  approximate 
evaluation  of  (68).  Note  that  ur,  urr,  etc.,  are  easy  to  find  at  r  ■  0. 


Consider  next  the  case  of  a  caustic.  A  caustic  will  have  curvature,  so  let 
the  caustic  be  a  circular  arc  (Fig.  6),  with  radius  R. 


(75) 


Again,  we  try  to  find  an  exact  solution  of  (63)  which,  in  the  ray  approximation, 
has  the  above  circular  arc  as  a  caustic.  We  try 

u(r,  6)  .  l  f(Y)  eilm(Y)"r  cos(8-Y»  dY  (76) 

0 

as  before.  Using  (67)  we  see  that  at  large  values  of  r  (0  <  0  <  T) ,  a  curve 
of  constant  phase  satisfies  r  +  m(0)  =  constant.  A  vector  tangential  to  this 
curve  has  (x,  y)  components  parallel  to  (-m*(0)  cos0  -  r  sin0, 

~m'  (0)sin0  +  r  cos0).  Since  a  vector  in  the  ray  direction  has  components 
(cosiji  ,  siniJO,  the  orthogonality  condition  becomes 

cost|»[“m'  (0)  cos0  -  r  sin0]  +  sinij)  [-m*(0)  sine  +  r  cos0]  =  0  (77) 

or 

-m'(0)  cos(0-ij/)  +  R(l-cost^)  =  0  (78) 

But  at  large  r,  0  -  ij»,  so  we  must  require 

m* (0)  =  R(l-cos40  (79) 


whence 


m(0)  =  R0  -  R  sin0 


(80) 


(the  constant  of  integration  is  absorbed  into  f(0)).  We  next  set 


m  -M  ei'/4 


(81) 


in  order  to  obtain  an  amplitude  at  a  large  reference  distance  A  of 


u(A,  6)  -  eik[R(6-cose)+A] 


(82) 


for  0  <  0  <r.  Then 
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(83) 


u(r,  0)  -/M  ei7r/4  /  eik[R(Y“sinY)  +  r  cos(0-Y)]  dy 
27T  0 

or  in  terms  of  p  and  \p, 

uCp.'iO  -Jf  eiir/4  /  eik [rY+r  sin(ij>-Y)  +  cosft>-Y)]  dy 

This  is  now  an  exact  solution  of  the  reduced  wave  equation;  at  large  r  the  ray 
approximation  yields  rays  having  the  desired  caustic. 


IVe  are  interested  in  u(p  ,iJ0  for  small  p,  i.e.,  close  to  a  caustic.  We  have  to 
break  up  the  calculation  into  two  regions.  In  the  first  of  these,  p  =  0(R); 
in  fact,  we  set  p  =  xR,  and  let  k  become  large  holding  x  fixed.  For  x  not  too 
large  and  for  \j)  not  too  close  to  0  or  T,  there  are  two  stationary  phase  points 
in  (84)  and  therefore 

u  -  ft  elkRt**T)  [1  ♦  leiRRO'0-*-2t)]  (8S) 

where  Yq  is  defined  by 


tan[Yo"^)/2]  *  x  (86) 

The  inverse  of  /T  implies  that  this  solution  breaks  down  near  x  *  0.  To 
handle  the  small  x  region,  we  use  a  boundary  layer  method.  Set 

u  =  a(x,  y)  eiWlO  (87) 

in  (63)  to  obtain 

lk(2cp  )opp  -  P-  op* 

+  ^2  “W  *  Cj  -  “1%  *  3  0  C88) 

If  we  discarded  all  terms  except  the  first  (for  large  k) ,  then  a  would  clearly 
be  singular  at  p  =  0.  We  therefore  set  p  =  k"s£|,  with  s  >  0;  a  substitution 
into  (88)  indicates  that  we  must  choose  s  =  1/3.  The  lowest  order  part  of  the 
boundary  layer  equation  thus  becomes 
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< 

» 

I 

) 


(89) 


il2“5  *  t”0  *  fl  ’  f?  “s  =  0 

whose  solution  has  been  given  in  terms  of  Airy  functions  by  Buchal  and  Keller 
(Reference  2)  as 


(1)  (2) 

u  =  av  J  +  c.a'-  J 


(90) 


„a>  .  expt-4)  Ai(2-2/3  ?  t-W) 


3R- 

«t2)  =  «p(-  Mj)  AU2-2/3  R-4/3  e4lr^3) 


(91) 


The  behavior  of  a  for  large  C  must  be  the  same  as  its  behavior  from  (85)  for 
small  x,  if  there  is  to  be  an  overlap  region.  After  some  algebra,  this 
requirement  yields: 


Cj  =  e-1*/12  25'6  4a  R-1/3  k1'6 


c2  =  i  e”'12  2S/6  4a  R*1'3  k1/6 


(92) 


In  particular,  the  value  of  u  on  the  caustic  itself  will  be 


6-1/6  a1/2  ,-1/2  r-1/3  ,1/6  r(l)  0iir/4 


(93) 
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3.  THE  ARCTIC  OCEAN  ENVIRONMENT  FOR 

/ 

UNDERWATER  ACOUSTIC  STUDIES 

/ 

The  principal  features  of  the  Arctic  Ocean  environment  of  concern  for  studies 
of  long  range  underwater  sound  propagation  are  bottom  topography  and  compo¬ 
sition,  ice  cover  roughness  distribution  and  sound  speed  profile  distribution. 
The  sound  speed  profile  is  in  turn  determined  by  the  local  temperature  and 
salinity  profiles.  Additional  areas  of  concern  for  acoustic  propagation 
studies  are  in  the  problems  of  propagation  losses  which  are  determined  by 
sea  water  absorption  and  scattering,  surface  (ice  cover)  reflection  and  absorp¬ 
tion  and  bottom  reflection  and  absorption. 


In  this  discussion  we  will  review  the  current  knowledge  of  these  topics  as 
available  from  tne  open  literature.  Extensive,  up  to  date  and  comprehensive 
information  on  acoustic  propagation  experiments  in  the  Arctic  Ocean  as  well 
as  on  topography  is  available  in  the  classified  literature.  The  two  references 
of  interest  here,  References  1  and  2,  deal  with  topography  and  transmission 
experiments  which  represent  the  current  knowledge  as  of  December  1972,  Ref¬ 
erence  3.  A  case  should  be  made  here  for  declassification  of  the  classified 
topographical  charts  to  make  more  comprehensive  and  uniform  data  available 
to  the  general  scientific  community.  \ 


3.1  BOTTOM  TOPOGRAPHY 

The  Arctic  Basin  may  be  divided  into  two  primary  topological  areas:  that 
comprised  of  continental  shelves  forming  about  two-thirds  of  the  total  area 
and  that  forming  the  central  basin.  The  continental  shelf  north  of  Eurasia 
is  the  widest  in  the  world  extending  up  to  850  km  width  in  the  Barents  Sea. 
This  shelf  is  cut  by  deop  re-entrants  north  of  Europe--the  best  known  are 
the  Svataya  Anna  and  Voronin  Troughs--and  a  marked  absence  of  troughs  north 
of  Asia.  North  of  Alaska  and  the  Chukchi  Sea  there  are  three  submarine 
canyons,  two  leading  into  the  Chukchi  Plain  and  one,  Barrow  Canyon,  leading 
into  the  Canada  Basin.  The  shelf  north  of  Canada  is  characterized  by  deep 
passages  among  the  islands  of  the  archipelago  leading  from  the  rather  rugged 
structural  highs  rising  from  the  floor  of  the  central  basin. 
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Hie  central  Arctic  is  divided  by  three  submarine  mountain  ranges  into  four 
basins.  These  mountain  ranges  run  roughly  parallel  to  oach  other  and  are 
aligned  to  150*  -  330°  longitudo.  The  central  range  is  the  Lomonosov 
Ridge  very  nearly  crossing  the  pole;  this  feature  is  very  narrow,  only  about 
40  km  wide  at  the  narrowest  point,  with  crest  depths  of  850  to  1200  meters. 

On  the  Alaskan  side  of  the  Ridge  lies  the  Marka-ov  Basin  (Wrangol  and  Fletcher 
abyssal  plains)  and  on  the  European  side,  the  Fram  Basin  (Pole  Abyssal  Plain). 
In  turn,  the  Fram  Basin  is  separated  from  the  Nautilus  Basin  (Barents  abyssal 
plain)  by  the  Nansen  Cordillera  (Mid-Oceanic  Ridge).  The  previously  noted 
Svataya  Anna  and  Voronin  Troughs  lead  into  the  Nautilus  Basin  between  tho 
surface  features  Severnaya  Zcmlya  and  Franz  Josef  Land.  Toward  the  Alaska- 
Canada  side  of  the  Lomonosov  Ridge  there  is  the  dominant  Alpha  Cordillera 
(known  as  the  Mendeleev  Ridge  in  USSR)  which  separates  the  Markarov  Basin 
from  the  extensive  Canada  Basin  (Beauford  Deep).  The  Alpha  Cordillera  is 
a  rather  broad  mountainous  feature  which  expands  to  join  the  Lomonosov  Ridge 
north  of  Greenland;  the  rugged  terrain  emanating  from  the  Canadian  Arctic 
Archipelago  can  be  thought  of  as  focusing  toward  the  Alpha  Cordillera  as 
one  moves  toward  the  USSR. 

The  greatest  density  of  bathymetric  data  is  available  for  the  sector  170°  to 
280°  longitude  bounded  by  the  US-Canada-Greenland  shoreline.  Data  have  been 
gathered  from  drifting  ice  island  stations  and  air  lifted  temporary  stations 
as  well  as  from  ship  and  submarine  operations.  In  a  recent  paper,  Wold  and 
Ostenso,  Reference  4,  present  a  detailed  bathymetric  chart  of  this  segment 
with  200  meter  contour  intervals  except  10-meter  contour  intervals  over 
Chukchi  shelf.  The  accuracy  of  the  bathymetry  data  resulting  from  the 
reported  survey  (from  1960  through  1969)  is  quoted  as  t  2%  with  location 
errors  of  i  3  km  for  the  ice  island  measurements  and  i  15  km  for  the  air 
lifted  station  measurements.  Hold  and  Ostenso  have  also  included 'hll  other 
soundings"  (Reference  4,  page  6254)  in  cogpiling  the  contour  map. 


*  Throughout  this  report  we  write  150"E  longitude  as  150*  and  30*W  as  330* 
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Fairly  detailed  bathymetry  is  also  available  for  the  Canadian  Arctic  from 
the  Canadian  charts,  Reference  5.  Since  the  sector  covered  includes  the 
sector  of  Wold  and  Ostenso,  these  charts  should  be  compared  for  consistency. 

A  comprehensive  physiographic  diagram  of  the  floor  of  the  Arctic  Ocean  was 
prepared  by  Dr.  Alan  Beal  for  the  U.  S.  Navy  Electronics  Laboratory  in  1966 
which  utilized  all  available  drift  station  soundings  in  addition  to  submarine 
soundings  through  1962.  Versions  of  this  data  are  available  in  Reference  6 
and  7.  However,  these  soundings  arc  concentrated  along  the  tracks  of  the 
drift  stations  and  submarines,  hence,  leaving  large  gaps  in  the  coverage 
where  the  information  must  be  extrapolated.  An  updating  of  this  information 
including  subsequent  submarine  soundings  has  been  in  progress  by  Beal  but, 
for  various  reasons,  is  now  at  a  standstill.  Reference  8. 

The  U.  S.  Naval  Oceanographic  Office  (NAVOCEANO)  has  bath'onetr’ c  charts 
covering  most  of  the  world's  oceans.  Since  NAVOCEANO  is  the  official  reposi¬ 
tory  of  all  sounding  data  obtained  by  ships  under  U.  S.  Navy  control  as  well 
as  collecting  data  from  other  Government  agencies  and  private  sources,  the 
NAVOCEANO  charts  are  also  a  good  source.  Specific  bottom  profile  accuracy 
can  be  improved  by  using  information  from  the  collection  sheets  (original 
sounding  echogram  entries)  compiled  Ly  the  Bathymetry  Division,  Hydrographic 
Surveys  Department,  NAVOCEANO,  which  are  used  in  constructing  the  topograph¬ 
ical  contour  maps. 

3.2  SOUND-SPEED  PROFILE  DATA 

Speed  of  sound  in  water  is  a  function  of  temperature,  salinity,  and  pressure. 
In  the  main  Arctic  basin  the  variation  of  temperature  and  salinity  is  such 
that  for  the  most  part  the  speed  of  sound  profile  is  an  increasing  function 
of  depth  from  the  surface  to  the  bottom.  Furthermore ,  as  noted  by  various 
investigators,  References  9  and  10,  the  acoustic  signal  propagation  charac¬ 
teristics  which  are  governed  by  the  sound-speed  profile  are  quite  stable 
both  in  time  and  geographic  location  indicating  that  the  temperature  and 
salinity  profiles  in  the  primary  Arctic  basin  do  not  vary  substantially. 


Variations  of  the  temperature  and  salinity  profiles  can,  however,  be  expected 
near  the  Bering  Strait  and  junctions  with  the  Greenland  and  Norwegian  Seas 
where  the  influence  of  Pacific  and  Atlantic  Ocean  waters  is  most  significant 
(References  11,  12,  and  13).  As  a  result  of  the  characteristic  sound-speed 
profiles,  there  exists  a  natural  half-channel  wave  guide  in  the  Arctic  with 
the  axis  at  the  undericc  surface  or  surface  of  open  water.  This  Arctic  wave¬ 
guide  is  the  low  frequency  analog  of  the  deep  sound  channel  or  SOFAR  channel 
of  the  non-polar  oceans  (low  frequency  because  the  underice  roughness  causes 
loss  of  high  frequency  signals  after  some  distance  of  propagation). 

The  prime  source  of  sound-speed  profile  data  in  the  oceans  is  the  National 
Oceanographic  Data  Center  (NODC) ,  Reference  14.  The  data  is  stored  on 
magnetic  tape  and  consists  of  temperature,  salinity  and  sound  speed  versus 
depth.  Sound  speed  is  calculated  by  Wilson's  equation,  Reference  15. 

The  NODC  data  file  can  be  used  to  generate  the  sound  speed  field  along  any 
desired  propagation  path  by  selecting  a  profile  data  zone  (determined  by  a 
strip  extending  to  each  side  of  the  propagation  path  to  some  distance  at 
which  profiles  arc  to  be  considered  satisfactory  for  inclusion  in  actual 
path  profile  determination).  A  search  program  operated  by  NODC  will  then 
search  the  magnetic  tape  and  print  out  those  sound  speed  profiles  taken  in 
the  month  or  season  of  interest  and  located  in  the  selected  data  zone. 

Data  in  the  file  is  keyed  by  the  Marsden  Square  Chart. 

Data  on  the  microstructure  of  the  sound-speed  profile,  as  determined  by 
temperature,  salinity  and  pressure  variations  is  not  readily  available. 
Current  work  under  the  AIDJEX  project  by  Coachman,  et  al.,  Reference  18,  is 
yielding  a  more  detailed  understanding  of  possible  layering  in  the  upper 
regions  of  the  Arctic  Ocean  water  mass.  However,  since  the  AIDJEX  stations 
cover,  at  best,  a  200  x  200  km  area  of  the  Ocean  (100  km  station  to  station 
distance  plus  maximum  drift  expectation),  the  data  may  not  be  of  immediate 
use  in  long  range  acoustic  signal  propagation  studies.  This  data  enters 
the  NODC  file  with  about  a  year's  delay.  From  the  point  of  view  of  the  low 
frequencies  of  interest  (~  20  Hz),  the  microstructure  will  not  have  signi¬ 
ficant  effects. 
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3.3  UNDER ICE  CHARACTERISTICS 


Since  one  boundary  of  the  sound  channel  is  formed  by  the  ice  canopy  (or,  on 
occasions,  the  Arctic  open  water  surface),  the  characteristics  of  the  under¬ 
ice  surface  must  be  known.  The  sound  waves  propagate  by  refraction  in  the 
deep  water  and  repeated  reflections  from  the  ice  canopy  so  the  primary  feature 
of  interest  is  the  underice  roughness.  The  dominant  roughness  features  are 
in  the  form  of  pressure  ridges  with  above  and  below  ice  buildup  due  to  the 
motion  and  resultant  deformation  of  the  Arctic  ice  pack.  Information  on  the 
characteristic  dimensions  of  the  above-ice  features  (termed  "sails")  and 
underice  features  ("keels")  has  been  derived  from  overflight  observations 
using  visual  and  laser  profilometer  techniques  and  submarine  sonar  soundings. 
Such  information  allows  a  compilation  of  sail  height  and  keel  depth  distri¬ 
butions  and  also,  to  some  extent,  the  distribution  of  spacings  between  the 
ridges.  The  significant  papers  on  this  data  are  in  References  19,  20,  and  21. 

Since  the  distribution  of  the  pressure  ridges  is,  to  current  knowledge, 
random,  a  statistical  description  of  the  underice  roughness  is  most  suited 
for  sound  propagation  studies.  Reference  20  derives  a  theoretical  distribution 
function  for  pressure  ridge  sail  heights  and  keel  depths  from  the  assumption 
of  randomness  and  also  presents  a  distribution  function  for  ridge  spacings 
based  on  spatial  randomness  assumption.  These  distrib  ition  functions  are 
correlated  with  observations  with  good  agreement.  From  the  distribution 
functions,  values  for  the  rms  ice  roughness  may  be  determined  to  be  used 
directly  in  sound  propagation  models.  Further  detail  on  the  underice  rough¬ 
ness  will  be  available  for  limited  regions  of  the  polar  pack  from  underice 
profilometer  experiments  currently  in  progress  in  conjunction  with  the  AIDJEX 
program  by  the  Applied  Physics  Laboratory,  University  of  Washington,  Seattle 
(Mr.  Robert  Francois).  This  data  has  not  yet  been  correlated  with  the  data 
and  theory  of  References  19  and  20. 

3.4  PROPAGATION  LOSS 

To  compute  the  sound  intensity  at  any  listening  location,  one  must  take  into 
account  the  propagation  losses  due  to  absorption  in  the  water  and  by  reflections 
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from  the  boundaries  (ice  cover  and  bottom)  of  the  sound  channel.  The  absorp¬ 
tion  loss  for  the  low  frequencies  of  interest  in  long  range  propagation 
studies  is  small;  however,  it  may  be  calculated  by  extrapolation  of  the  Marsh 
and  Schulkin  formula  of  Reference  22. 

Thorp,  Reference  23,  has  advocated  an  alternate  (simpler)  formula  for  attenu¬ 
ation  at  low  frequency;  however,  in  a  review  paper  (Reference  24)  published 
in  March  1972,  Berman  and  Guthrie  make  the  observation  that 

"...  we  find  that  when  such  a  critical  analysis  (of  experimental 
evidence)  is  made,  little  or  no  valid  evidence  exists  to  point 
to  the  existence  of  a  low  frequency  absorption  above  that  which 
may  be  extrapolated  from  the  Marsh  and  Schulkin  formula." 

Furthermore,  in  view  of  the  uncertainties  of  the  boundary  characteristics 
and  their  contribution  to  propagation  loss,  the  differences  in  the  loss  coef¬ 
ficients  are  negligible. 

Exceedingly  little  data  exists  on  bottom  reflection  losses  for  the  calculation 
of  which  bottom  density  and  sound  speed  are  required.  In,  the  review  paper, 
Berman  and  Guthrie,  speaking  about  ocean  propagation  in  general,  state  that 

"Scant  data  exists  on  bottom  reflection  loss  to  serve  as  a  guide 
for  even  representative  models  for  the  low  frequency,  low  grazing 
angle  reflections  that  can  be  expected  to  dominate  in  long  range 
propagation  by  bottom  bounce.  It  can  be  assumed  that  a  ray  inci¬ 
dent  on  the  bottom  will  undergo  a  fractional  attenuation  that  is 
a  function  of  wavelength,  grazing  angle,  surface  roughness  and 
the  geological  structure  of  the  bottom." 

Information  on  the  geological  structure  of  the  bottom  could  become  more  avail¬ 
able  with  the  explorations  by  commercial  interests  for  mineral  resources  in 
the  Arctic  spurred  by  the  oil  discoveries.  However,  much  of  this  information 
will  never  be  available  to  the  scientific  community  unless  specific  coordi¬ 
nation  with  commercial  exploration  efforts  is  made  by  interested  Government 
agencies. 
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The  most  up  to  date  and  comprehensive  results  on  experimental  studies  of 
propagation  losses  in  the  Arctic  are  those  given  in  the  classified  report 
by  Buck,  et  al.,  Reference  2.  Considerable  relative  loss  data  have  been 
gathered  by  Kutschale  and  co-workers  at  Lamont- Doherty  Laboratories , 

Columbia  University,  New  York,  which  are  useful  in  comparative  acoustic 
propagation  studies. 

3.5  RECOMMENDATIONS  FOR  ACQUISITION  OF  REQUIRED  DATA 

The  review  of  available  (unclassified)  literature  indicates  that  substantially 
greater  effort  is  necessary  to  acquire  the  necessary  density  and  quality  of 
environmental  data  of  importance  to  long  range  sound  propagation  in  the  Arctic. 
Most  critical  shortcomings  exist  in  the  definition  of  bottom  topology  and  geo¬ 
logical  structure  as  well  as  underice  characteristics.  The  density  of  bathy¬ 
metric  data  as  well  as  the  quality,  can  be  greatly  increased  at  reasonable 
cost  by  use  of  a  recently  developed  airlift  through  the  ice  bathymetric  system. 
Edo  Western  Corporation,  Reference  25,  has  tested  a  prototype  transducer  which 
can  be  placed  on  the  snow  or  ice  and  operated  by  remote  control  from  inside  a 
helicopter  or  ski  plane.  The  battery  powered  system  generates  a  high  power, 
low  frequency  signal  capable  of  penetrating  ice  of  up  to  12  feet  in  thickness. 
Depth  readings  of  1400  feet  were  encountered  in  the  prototype  test  and  typical 
time  per  data  point  was  on  the  order  of  one  minute.  Development  of  rapid 
salinity  and  temperature  sampling  procedures  would  allow  thorough  mapping  of 
the  medium.  Detailed  low  frequency  propagation  loss  studies  should  be  carried 
out  in  conjunction  with  the  oceanographic  data  gathering  along  specified  known 
paths  using  harmonic  sources.  Use  of  underice  profilometers  along  the  selected 
propagation  paths  would  define  more  precisely  the  effects  of  underice  roughness 
on  the  signals. 
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4.  APPLICATION  OF  RAY  ACOUSTICS  IN  THE  ARCTIC 

The  factors  that  affect  underwater  sound  propagation  are  the  nature  of  the 
sound  speed  depth  profile,  bottom  and  surface  reflections  and  losses  due  to 
spreading  and  absorption.  The  sound  speed  depth  profile  is  determined  by 
the  vertical  variation  of  salinity,  temperature  and  pressure.  Bottom  reflec¬ 
tion  losses  are  dependent  on  the  nature  of  the  bottom  sediment ,  surface 
roughness  and  topography.  Surface  reflection  losses  are  functions  of  ray 
^razing  angle  and  the  roughness  of  the  ice  cover.  Absorption  losses  are 
due  to  the  internal  generation  of  heat  as  sound  propagates.  All  the  ref¬ 
lection  and  absorption  losses  are  very  much  dependent  on  the  frequency  of 
propagation. 

In  the  Arctic  it  is  known  that  below  the  upper  layer  of  several  hundred 
feet  the  sound  speed  gradient  with  respect  to  the  depth  is  fairly  uniform 
and  moderate.  In  addition,  measurement  data  tend  to  indicate  that  spatial 
variation  of  such  vertical  sound  speed  profiles  are  small,  especially  in  the 
central  Arctic  region.  Under  these  conditions  ray  theory  affords  a  very 
convenient  and  effective  technique  for  calculating  the  acoustic  intensities, 
ruling  out  those  cases  in  which  extreme  low  frequencies  are  considered 
(Section  2.2). 

A  typical  ray  diagram  and  sound  speed  profile  using  recent  temperature  and 
salinity  data  at  Base  Camp  Station  (72-0, 20N  latitude  and  211-24.8  longitude) 
are  plotted  and  presented  in  Section  6.  The  units  of  depth  and  range  are 
in  feet  and  the  speed  of  sound  in  feet  per  second.  The  source  is  located 
1000  feet  below  the  surface  and  the  initial  grazing  angles  of  the  rays 
have  one  degree  increment. 

Such  ray  diagrams  form  the  basis  of  calculating  the  acoustic  intensities 
by  ray  theory.  A  program  has  been  developed  that  is  designed  specifically 
for  such  sound  intensity  calculations  at  a  given  number  of  receiver  stations 
corresponding  to  a  given  number  of  sources.  The  locations  of  the  receiver 
stations  and  sources  are  specified  in  terms  of  latitude,  longitude,  and 
depth.  Input  parameters  to  the  program  Are  the  frequency  of  propagation, 
typical  salinity  and  temperature  profiles,  bottom  topography  and  the  RMS 
ice  roughness.  ,c 


4.1  SOUND  SPEED  PROFILE  CALCULATION 


The  sound  speed  is  a  function  of  the  local  salinity,  temperature  and 

pressure.  In  general  it  is  calculated  by  the  Wilson’s  formula 
(Reference  1) : 

v  =  1449.14  +  vT  +  vp  +  vs  +  Vgjp 

where 

v  =  4.5721T  -  4.4532X10"2T2  -  2. 6045X10_4T3  +  7.9851X10‘6T4 
vp  a  1.60272X10'1P  +  1.0268X10"5P2  +  3. 5216X10‘9P3  -  3. 3603X10“12P4 
vg  =  1 . 39799 (S-35)  +  1 .69202X10"3(S-35) 2 
VSTP=  CS-35) C-l. 1244X10-2T  +  7. 7711X1(T7T2  +  7.7016X10_5P 
-  1 .2943X10“ 7P2  +  3. 1580X10“8PT  +  1 .5790X10“9PT2)  + 
P(-1.8607X10-4T  +  7.4812X10_6T2  +  4.5283X10-8T3)  + 

P2(-2.5294X10“7T  +  1.8563X10-^2)  +  P3(-l.  9646X10"  10T) 

T  is  the  temperature  in  degrees  Centigrade, 

P  is  the  pressure  in  kg/cm2  (absolute) , 

S  is  the  salinity  in  °/oo,  and 
v  is  the  sound  speed  in  m/sec. 

The  pressure  at  any  depth  z  may  be  calculated  by  Leroy's  formula 
(Reference  2) : 

P  =  1.04  +  0.102506(1  +  0.00528  sin2<#>) z  +  2.524X10“7z2  (2) 


where 

z  is  the  depth  in  meters  and 
<{>  is  the  latitude  in  degrees. 

Although  the  Wilson's  formula  is  generally  used  for  sound  speed  calculations, 
a  shorter  and  simpler  version,  claimed  to  fit  applicable  Wilson's  data  better 
is  due  to  Leroy  (Reference  3) : 


v0  =  1493.0  +  3(T-10)  -  6X10'3(T-10)2  -  4X10"2(T-18) 2 

+  1.2 CS- 35)  -  10"2(T-18) (S-35)  +  z/61 

va  =  lO’1?2  +  2X10"4S2(T-18)2  +10-1?4*/90 

v  =  2.6X10“4(T-5)  (T-25)T 
b 

v£  =  -10"V(C-4)(?-8) 

vd  =  1.5X10"3(S-35)(l-£)  +  3X10_6T2 (T-30) (S- 35) 

T  is  the  temperature  in  degrees  Centigrade, 

S  is  the  salinity  in  °/oo, 

<t>  is  the  latitude  in  degrees , 
z  is  the  depth  in  meters,  and 
Z,  =  z/1000,  i.e.,  depth  in  km. 

Leroy  claims  that  not  all  data  used  in  establishing  Wilson's  formula  are 
realistic  in  view  of  the  actual  conditions  to  be  expected  in  the  natural 
sea  waters  of  the  world  and  uses  what  he  considers  a  more  representative 
data  base  to  establish  his  formula. 

To  find  the  sound  speed  at  a  particular  depth  corresponding  to  the  given 
data  of  salinity  and  temperature  at  various  depths,  one  first  calculates 
the  salinity  and  temperature  at  the  desired  depth  by  spline  interpolation . 
This  ensures  "smooth"  salinity-depth  and  temperature-depth  profiles 
through  the  given  data  points  in  the  sense  not  only  the  first  but  also 
the  second  derivatives  are  continuous  from  interval  to  interval.  From 
the  interpolated  values  of  salinity  and  temperature  the  sound  speed  is 
then  calculated  by  either  the  Wilson's  or  Leroy's  formula. 

In  order  to  compare  Wilson's  and  Leroy's  formulas  sound  speed  has  been 
calculated  at  different  depths  (using  spline  interpolation  of  salinity 
and  temperature  as  described  above)  up  to  12000  feet  in  depth  (Base 
Camp  Station  data  of  salinity  and  temperature) .  The  results  indicate 
that  differences  in  sound  speed  calculated  by  both  methods  are  insigni¬ 
ficant,  being  of  the  order  of  1/10  of  one  percent  at  larger  depth.  It 
thus  appears  that  Leroy's  formula  can  indeed  be  used  to  advantage  because 
of  its  relative  simplicity. 
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4.2  BOTTOM  TOPOGRAPHY 


In  order  to  give  a  fairly  adequate  description  of  the  topography  of  the 
Arctic  Ocean  bottom,  4800  stored  words  are  assigned  for  this  purpose. 

This  covers  the  rectangular  area  of  which  the  center  is  the  North  Pole: 
the  midpoints  of  the  four  sides  have  the  following  latitudes  and  longi¬ 
tudes:  70°N  0°,  75°N  90°,  70°N  180°  and  75°N  270°.  This  rectangle  is 
partitioned  into  4800  squares  by  dividing  the  horizontal  sides  into  80 
equal  intervals  and  the  vertical  sides  into  60  intervals.  Each  such 
square  thus  has  a  side  equal  to  about  34.55  miles,  an$  each  stored 
word  represents  the  average  depth  of  the  ocean  bottom  for  this  area. 

Although  a  number  of  sources  are  available  for  the  bathymetric  data  of 
the  Arctic  Ocean,  it  is  finally  decided  to  make  use  of  the  Arctic  Bathy¬ 
metric  maps  as  compiled  by  the  Canadian  Hydrographic  Service,  1966-1967, 
which  furnish  fairly  detailed  information.  However,  these  maps  cover 
only  the  two  quadrants  of  the  Arctic:  latitude  72°N  longitude  180° -270° 
and  latitude  72°N  longitude  270°-360°;  viz.,  these  areas  on  the  side  of 
Alaska,  Canada  and  Greenland.  Fortunately  much  information  can  also 
be  derived  from  the  North  Polar  Chart  compiled  by  the  Hydrographic 
Department  of  the  British  Navy,  1969,  which  supplies  the  detailed  bathy¬ 
metric  data  on  the  Russian  and  European  side. 

A  bathymetric  computer  plot  showing  the  approximate  depth  of  the  Arctic 
Ocean  bottom  at  4800  locations  is  given  in  Section  6.  As  noted  above, 

each  letter  or  numeral  represents  the  average  depth  of  an  area  of  about 
34.55  mile  square. 

In  order  to  give  a  visual  description  of  the  bottom  topography,  a  number 
of  Arctic  Ocean  contour  plots  are  presented  in  Section  6,  based  on  the 
stored  depth  data  and  plotted  by  CALCOMP.  These  plots  give  the  various 
depth  contours  of  100,  300,  500,  700,  900  ft.;  1000,  2000,  3000,  4000, 
5000  ft.;  6000,  7000,  8000,  9000,  10000  ft.  and  finally  11000,  12000, 
13000,  14000,  15000  ft.  In  each  of  these  plots  the  shore  lines  are  also 
included.  Although  not  explicitly  indicated  in  these  plots,  the  Lomono¬ 
sov  Ridge,  the  Canadian  Abyssal  Plain,  the  Amundsen  Basin,  the  Fletcher 

Plain,  etc.,  are  distinctly  visible. 
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4.3  RMS  ICE  ROUGHNESS 


While  it  is  possible  to  determine  the  RMS  ice  roughness  from  statistical 
considerations  of  the  ice  ridge  spacings,  sail  heights  and  keel  depths 
(Section  3.4),  the  acoustic  experiments  of  Mellen  and  Marsh  (Reference  4) 
and  Kutschale  (Reference  5)  directly  relate  the  ridge  index  (number  of 
ridge  counts  per  30  nautical  miles)  with  the  RMS  ice  roughness.  These 
experiments  indicate  the  following  relationship: 


Ridge  Index 


RMS  Ice  Roughness 


300  -  600 
600  -  900 
900  -  1200 


2. 5  meters 

3.0 

4.0 


The  results  of  computation  in  this  report  are  based  on  tho  ridge  index 
data  in  the  spring  season  supplied  by  Whittmann  and  Schule  (Reference  6) . 

A  plot  of  RMS  ice  roughness  utilizing  the  above  relationship  is  presented 
in  Section  6. 


4.4  ICE  REFLECTION  LOSS 

It  is  reasonable  to  assume  that  for  grazing  angles  less  than  30°,  little 
energy  is  lost  through  penetration  into  ice.  For  long  range  propagation 
only  those  rays  that  have  perfect  reflection  at  the  underice  surface 
have  significant  contribution  to  the  acoustic  intensities.  Marsh  (Re¬ 
ference  7)  has  shown  that  the  loss  per  such  reflection  can  be  calculated 
by  the  following  equation: 

a  =  -10  log1(J(l  -  2.57X10'4f3/2h8/5  sin6  )  (4) 

where  a  is  the  loss  in  db  per  reflection, 
f  is  the  frequency  in  hz., 
h  is  the  RMS  ice  roughness  in  meters  and 
0  is  the  grazing  angle  at  incidence  of  limiting  ray. 

As  can  be  seen  from  (4),  the  reflection  loss  increases  rapidly  with 
frequency. 
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4.5  A3S0RPTI0N  LOSS 

This  is  the  loss  due  to  the  conversion  of  acoustic  energy  into  heat 
when  sound  propagates  in  a  medium.  It  can  be  calculated  according  to 
the  formula  suggested  by  Marsh  and  Schulkin  (Reference  8) : 

/  ASf?f2  Bf2  ^  A 

0  S(lf7  f 2  +  ~  6.54X10"4P)(7.943X10-3)  (* 

where  0  is  the  absorption  loss  coefficient  in  db  per  kyd, 

A  =  2.34X10  6,  constant  for  ionic  relaxation  in  sea  water. 

-6  1 
B  =  3.38X10  ,  constant  associated  with  pure  water  viscosity 

mechanism, 

S  is  the  salinity  in  °/oo, 

P  is  the  pressure  in  kg/cm2, 

T  is  the  temperature  in  degree  Centigrade, 
f  is  the  acoustic  frequency  in  kcps,  and 
fT  =  2. 19X10^  1520/(T+273) )  ^  tjje  temperature  dependent 
relaxation  frequency  in  kcps  at  atmospheric  pressure. 


For  the  low  frequency  under  consideration,  viz.,  20  hz. ,  the  above 
formula  has  to  be  extrapolated.  An  alternate  formula  for  calculating 
the  absorption  loss  coefficient  is  due  to  Thorp  (Reference  9) : 

O.lf2  40 f 2 

1  +  f2  4100  +  f2 

Formula  (6)  is  obtained  by  fitting  curves  to  actual  measured  data  of 
low  frequencies.  For  the  same  f  (6)  gives  much  higher  values  of  0 
than  those  extrapolated  by  using  (5).  For  example,  for  f  =  ,020(20hz.), 

5  =  34,  T  =  0°c,  0  a  4.3886X10'S  db/kyd  as  calculated  by  (6)  and 

6  =  4.29698X10  11  db/kyd  by  (5).  But  even  using  the  larger  value  of  0 
the  transmission  loss  is  only  0.048  db  for  sound  travel  of  1000  km. 

The  absorption  loss  is  therefore  insignificant  for  the  low  frequency 
under  consideration. 


Formula  (5)  is  used  in  the  program  and  throughout  all  computations. 

However,  in  view  of  the  large  discrepancies  between  the  results  calculated 

by  (5)  and  (6) ,  further  wo1"’'  remains  to  done  in  this  area  when  transmission 
of  higher  frequencies  is  considered. 
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4.6  BOTTOM  REFLECTION  LOSS 


A  detailed  computation  of  the  bottom  reflection  loss  has  not  been  carried 
out  for  the  following  reason.;: 

(a)  The  bottom  reflection  loss  is  a  function  of  the  grazing  angle, 
the  frequency  of  propagation,  surface  roughness  and  the  nature 
of  the  bottom  sediment.  Data  on  the  latter  are  very  scarce. 

(b)  For  long  range  propagation  those  rays  that  "bounce"  many  times 
on  the  bottom  before  reaching  the  receiver  will  suffer  enough 
attenuation  and  their  contribution  to  the  sound  intensity  is 
insignificant  compared  with  those  rays  that  reach  the  receiver 
without  touching  the  bottom. 

(c)  To  trace  each  ray  from  the  source  to  the  receiver  over  the 
bottom  topography,  which  is,  at  best,  some  approximation, 
requires  enormous  amount  of  computing  time. 

It  therefore  stands  to  reason,  at  least  for  a  preliminary  study,  to  con¬ 
sider  only  those  rays  which  travel  in  a  uniform  channel  bounded  by  the 
surface  and  an  artificial  flat  bottom  of  some  effective  depth.  Those 

rays  which  penetrate  deeper  than  the  effective  depth  are  excluded  from 
consideration. 


Fig.  1 
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*0  determine  the  effective  depth  dg  from  source  S  to  receiver  R  one  finds 
the  successive  depths  dj,  d2,  dn  from  the  stored  depth  data  by 
tracing  a  ray  from  S  to  R.  Then  dg  can  be  calculated  in  one  of  the 
following  ways: 


de  »  min  (d1(  d2,  — ,  dn) 


(7A) 

(7B) 


As  shown  in  the  figure  dQ  occurs  at  X  based  on  (7A)  but  some  rays  can  still 
travel  from  S  to  R  without  touching  the  bottom.  Eq.  (7B)  is  probably  more 
realistic  because  it  is  the  RMS  of  the  reciprocals  of  the  depths  and  gives 
emphasis  to  the  minimum  depth.  This  is  the  formula  used  in  all  the  compu¬ 
tations  to  determine  an  arbitrary  bottom  of  the  sound  propagation  channel. 


Since  the  maximum  grazing  angle  at  the  ice  cover  is  limited  to  30°  (Section 
4.4),  and  since,  by  assumption,  those  rays  penetrate  below  the  effective 
depth  do  not  contribute  to  the  acoustic  intensity  at  the  receiver,  all  the 
rays  to  be  traced  are  enclosed  between  the  two  extreme  rays  making  the  apex 
angle  A0  at  source  S  as  shown. 
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4.7  STATISTICAL  SUMMATION  OF  RAY  INTENSITIES 


The  intensity  I  of  each  ray  at  range  L  from  the  source  is  given  by: 


I  =  -P 


cos  0 


3L 

30, 


sin  0 


(8) 


where  I  is  the  intensity  in  energy  units  per  second  per  unit  area, 

P  is  the  source  strength  in  energy  units  per  second, 

0O  is  the  initial  grazing  angle  of  the  ray  at  the  source,  and 
0  is  the  grazing  angle  of  the  ray  at  the  point  under  consideration. 
Thus  to  determine  I  for  each  ray  an  "adjacent"  ray  at  grazing  angle 
0o+A0q  has  to  be  traced  simultaneously  in  order  to  evaluate  3ty390. 


To  calculate  the  acoustic  intensity  at  receiver  R  at  given  range  L  from 
the  source  S,  one  should  theoretically  find  all  the  rays  from  S  that  pass 
through  R  and  sum  up  the  intensities  of  the  individual  rays,  taking  into 
account  the  phase  changes,  etc.  In  view  of  the  fact  that  the  path  of  a 
ray  is  greatly  dependent  at  any  point  of  its  travel  on  the  speed  of  sound 
which  is  only  approximately  known,  any  slight  variation  in  the  speed  of 
sound  can  cause  the  ray  to  deviate  considerably  from  reaching  the  receiver 
R.  Such  a  procedure  is  therefore  only  of  academic  interest,  and  the 
following  one  is  used  instead  which  is  also  computationally  more  efficient. 


L 


We  trace  all  the  rays  from  S  enclosed  between  the  two  extreme  rays 
(Section  4.  6)  between  the  upper  and  lower  boundaries  with  successive 
increment  of  6  in  the  initial  grazing  angles.  Thus  if  N  rays  are 
traced,  6  =  A0/N. 

Now  at  distance  L  from  S,  these  rays  will  terminate  at  different  depths 
di  (i  =  1,  2,  ,N) ,  which  are  in  general  different  from  the  depth  of 

the  receiver,  d^.  However,  because  the  speed  of  sound  is  only  approxi¬ 
mate,  each  ray  has  a  chance  of  reaching  R,  and  we  assume  that  the  pro¬ 
bability  of  the  ith  ray  to  terminate  at  the  receiver  depends  on  the 
distance  |dR-dJ  non-dimensionalized  by  dfi,  the  effective  depth,  i.e., 

I dR-di I /de  and  is  given  by  the  normal  distribution: 


1 

— _ —  e 
a/ 2tt 


(dn-di) 

2d'i  a* 
e 


(9) 


where  o  is  the  standard  deviation.  If  only  a  small  number  of  rays  are 
traced,  obviously  the  bell  shaped  normal  distribution  should  be  flatter 
to  include  more  rays  at  R,  so  it  is  reasonable  to  assume  that  a  depends 
on  the  number  of  rays  traced:  a  =K|/vfr,  k, being  an  input  parameter  to  the 

program  (Reference  10).  The  average,  or  expected  value  of  the  intensity 
at  R  is  therefore 


_ 1_ 

o/2tt 


Z  I  e 
i 


tyi>2 


(10) 


where  Ii  is  the  intensity  of  the  ith  ray  at  range  L  from  S.  Note  that 
incoherent  summation  of  the  rays  (phase  neglected)  ±s  used  here. 


Since  most  measurement  data  are  expressed  in  transmission  loss  (db) ,  for 
purposes  of  comparison  the  corresponding  transmission  losses  are  computed 
according  to  the  formula:  10  log10(Io/I),  where  IQ  is  the  intensity  at 
1  yard  from  the  source. 
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4.8  CORRECTION  FOR  CAUSTICS 

Abnormal  high  intensities  of  some  of  the  rays  may  occur  if  for  a  given 
range  from  a  source  a  position  is  reached  which  may  be  close  to  a  caustic 
zone.  Such  a  situation  is  likely  to  occur  especially  when  a  large  number 
of  rays  are  traced  and  will  result  in  undue  increase  in  the  acoustic  in- 
tensity  when  calculated  by  (10). 

in  order  to  correct  this,  the  intensities  of  all  rays  traced  are  made  to 
go  through  a  filter  to  eliminate  those  values  which  markedly  exceed  the 
average  value.  The  filtering  device  then  detects  these  abnormal  large 
values  of  f.  •  |l.-IA|/IA,  «h«e  IA  is  the  arithmetic  mean  of  It, 

I  =  EI./N  N  being  the  total  number  of  rays  traced.  Normally  fj  should 
be  less'than  1;  if  ^  >  *2.  where  *2  is  an  input  parameter  to  the  program 
(say  <2  -  2),  Ij  will  be  considered  to  be  close  to  a  caustic  and  is 
simply  rejected. 
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4.9  LLOYD  MIRROR  CORRECTION 


For  shallow  sources  or  receivers  the  phase  difference  of  the  rays  that 
arrive  directly  and  by  surface  reflection  will  result  in  a  modulation 
factor  of  the  intensity  given  by  4  sin2(  (27id/X)  sin9)  ,  where  d  can  be 
either  the  receiver  depth  or  the  source  depth,  X  the  wave  length  and  0 
is  the  ray  arriving  grazing  angle.  This  factor  is  applied  to  each 
individual  ray  traced,  and  is  not  taken  into  consideration  (i.e.,  the 
factor  is  equal  to  1)  if  d  exceeds  a  certain  depth,  usually  k-X,  where 
<~<1  and  is  an  input  parameter  to  the  program. 

4-10  ARCTIC  PARTITIONING  AND  SOURCE  AND  RECEIVER  LOCATIONS 

The  principal  mountain  ranges  in  Central  Arctic  are  the  Lomonosov  Ridge 
and  the  Arctic  Mid-Ocean  Ridge,  which  run  approximately  in  the  direction 
of  longitude  150°-330°.  Between  the  ridges  are  the  Amundsen  Basin  and 
the  Pole  Abyssal  Plain.  On  the  European  side  of  the  Arctic  Mid-Ocean 
Ridge  is  the  Nansen  Basin  and  the  Barents  Plain;  on  the  Alaska  side  of 
Lomonosov  Ridge  is  the  Markarov  Basin  and  the  Alpha  Cordillera. 

The  sources  and  receivers  are  divided  into  five  groups.  The  basis  of 
the  partitioning  is  that  those  receiver  stations  of  one  group  can  monitor 
sources  of  the  same  group  much  better  than  those  sources  belonging  to  a 
different  group.  This  is  because  sources  and  receiver  stations  of  dif¬ 
ferent  groups  are  usually  separated  by  larger  distances,  and  possibly  by 
mountain  ridges  which  partially  obstruct  ray  propagation.  Any  movements 
in  a  deep  channel  can  best  be  monitored  by  receivers  in  that  channel. 

Group  I  consists  of  seven  sources  and  five  receivers  which  are  located 
in  the  sector  with  latitude  less  than  85°  and  longitude  from  0°  to  90°. 
This  covers  the  shallow  region  of  the  Barents  Sea,  the  Svyatata  Anna 
Trough,  part  of  the  Nansen  Basin  and  extends  to  the  side  of  Greenland. 

Group  II  consists  of  eight  sources  and  six  receivers  covering  the  Arctic 
Mid-Ocean  Ridge,  the  Amundsen  Basin  and  the  Lomonosov  Ridge.  In  general 
this  is  the  area  of  deep  waters. 


Group  III  consists  of  five  sources  and  three  receivers  covering  the  shallow 

region  of  Laptev  Sea  and  the  deep  Makarov  Basin  and  the  Fletcher  Plain.  It  is 
a  narrow  strip  which  runs  almost  parallel  to  the  Lomonosov  Ridge. 

Group  IV  consists  of  three  sources  and  three  receivers.  The  sources  are 

scattered  in  the  Alpha  Cordillera  and  the  receivers  are  close  to  Queen 
Elizabeth  Islands. 

Group  V  consists  of  four  sources  and  three  receivers  and  covers  the  East 
Siberian  Sea,  Chukchi  Plateau  and  the  Canada  Abyssal  Plain.  The  receivers 
are  close  to  the  shores  of  Alaska  with  receiver  number  1  very  near  the 
Naval  Arctic  Research  Laboratory  at  Point  Barrow. 

The  latitude  and  longitude  of  the  sources  and  receivers  of  the  various 
groups  are  given  in  tabular  form  in  Section  6.  Also  shown  are  the  sources 
and  receivers  in  the  Arctic  Ocean  Contour  Plots,  page  67  -  page  76. 

4.11  COMPARISON  OF  COMPUTED  RESULTS  WITH  FIELD  MEASUREMENTS 

For  each  partitioned  group  the  acoustic  intensity  and  transmission  loss 
are  computed  for  each  receiver  corresponding  to  each  sound  source.  The 
receivers  are  placed  at  depths  of  100,  200,  400,  600,  800,  1000,  1500, 

2000,  3000  and  4000  feet;  each  source  is  located  at  depths  of  20,  210, 

410  and  610  feet.  The  detailed  computed  results  are  presented  in 
Part  II  of  this  report.  The  range  is  the  horizontal  distance  between 
the  source  and  the  receiver  and  the  bottom  level  is  the  effective  depth 
calculated  according  to  Eq.  (7B) ,  Section  4.6,  between  the  source  and 
the  receiver,  both  expressed  in  feet.  The  frequency  of  propagation  is 
assumed  to  be  20  hz. 

.xtensive  measurements  of  sound  transmission  loss  in  the  Arctic  have  been 
carried  out  and  reported  by  Marsh  and  Mellen  (Reference  4)  and  also  by 
Buck  and  Greene  (Reference  11) .  These  measured  data  of  transmission  loss 
against  range  are  then  combined  and  smooth  curves  are  constructed  which 
bestfit  them  at  various  frequencies.  The  smooth  curves  express  the 
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relation  between  the  transmission  loss  with  range  based  on  considerations 
of  spreading  and  ice-water  reflection  loss  (Reference  12).  It  is  under¬ 
stood  that  in  the  experimental  measurements,  the  receivers  were  at  about 
60  meters  deep  and  the  sources  were  at  various  depths  up  to  250  meters. 

For  receiver  depth  of  200  feet  the  computed  transmission  loss  of  all  the 
cases  considered  is  plotted  against  range  irrespective  of  the  bottom  level. 
These  plots  are  given  in  Section  6,  page  77  -  page  88.  For  each  of  these 
plots  the  sources  are  located  at  one  of  the  depths:  20,  210,  410  or  610 
feet.  Thus  the  source  and  receiver  depths  correspond  roughly  to  the 
conditions  under  which  the  experimental  data  were  taken.  For  purposes 
of  comparison  the  best-fit  curve  given  by  Buck  (Reference  12)  for  20  hz. 
is  also  included  in  each  plot. 

It  is  interesting  to  note  that  for  shallow  sources  (20  feet  below  the 
surface)  the  calculated  results  all  fall  below  the  best-fit  curve  by 
Buck.  The  associated  large  transmission  loss  could  be  attributed  to 
the  Lloyd  mirror  effect.  For  sources  210  feet  below  the  surface  the 
calculated  results  are  above  the  best-fit  curve  while  for  sources  410 
and  610  feet  below  the  surface,  the  calculated  results  again  fall  below 
this  curve.  Thus,  on  the  average,  the  computed  results  agree  quite 
well  with  experimental  data. 

4.12  COMPILER  PROGRAM 

The  ray  tracing  program  is  designed  for  the  following  purposes: 

(a)  From  the  given  location  of  sources  and  receivers  the  program 
first  finds  the  depth  contours  (depth  values  at  regular  inter¬ 
vals)  from  each  source  to  each  receiver,  then  calculates  the 
effective  depth  according  to  (7B) ,  Section  4.6. 

(b)  The  program  then  calculates  the  acoustic  intensity,  and  the 
transmission  loss  from  each  source  to  each  of  the  receivers 
corresponding  to  the  specified  depths. 

The  program  actually  consists  of  two  parts.  From  the  input  data  the 
first  part  performs  (a) ;  the  calculated  depth  data  are  combined  with  the 

input  data  and  fed  into  the  second  part  to  perform  (b) . 
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The  execution  time  depends  greatly  on  the  number  of  rays  traced  in  each  case. 
For  tracing  on  the  average  of  90  rays  from  one  source  to  one  receiver,  it  takes 
about  160  seconds  of  CP  time  on  the  CDC  6600,  or  about  30  seconds  of  CP  time 
on  the  CDC  7600.  The  core  required  is  65K  with  the  depth  data  written  on 
tape. 
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5.  SUBCONTRACT  TO  LAMONT- DOHERTY 


The  subcontract  to  Lamont -Doherty  calls  for  the  development  of  appropriate 
techniques  to  calculate  the  acoustic  intensities  in  the  Arctic  Ocean,  and 
comparisons  of  computed  results  by  ray  theory,  normal  mode  theory,  or  any 
newly  developed  methods  against  experimental  data.  The  principal  effort 
is  thus  to  search  for  a  rapid,  accurate  method  of  computing  propagation 
loss  as  a  function  of  range  in  the  ice  covered  Arctic  Ocean.  Important 
input  parameters  to  the  propagation  model  are  ice  roughness,  bottom  topo¬ 
graphy,  and  the  velocity  structure  as  a  function  of  depth  in  the  ice,  water, 
and  bottom  sediments.  Computational  speed  is  of  the  utmost  importance  to 
evaluate  variations  of  these  parameters  on  propagation  loss  as  a  function 
of  range  since  a  number  of  models  must  be  considered. 

Such  a  rapid  computational  method  is  the  Fast  Field  Program  (FFP)  technique 
• discovered  by  Marsh  (1967) .  The  FFP  was  implemented  for  computation  on  a 
digital  computer  by  Di  Napoli  (1971)  for  propagation  in  an  all  liquid  wave 
guide.  In  the  Arctic  Ocean  solid  layers  as  well  as  liquid  layers  must  be 
considered,  and  a  convenient  starting  point  for  the  FFP  is  the  integral 
solution  derived  by  matrix  methods  by  Kutschale  (1970,  1972)  for  propagation 
from  point  harmonic  sources  in  a  multilayered  liquid-solid  half  space. 

In  the  FFP  technique  the  integral  solution  is  evaluated  rapidly  as  a  func¬ 
tion  of  range  by  numerical  integration  employing  the  Fast  Fourier  Transform 
(FFT).  Singularities  in  the  integrand  corresponding  to  the  normal  mode 
poles  are  removed  from  the  axis  of  integration  by  including  an  attenuation 
factor  in  each  layer.  The  attenuation  in  each  layer  can  represent  the 

effects  on  propagation  loss  as  a  function  of  range  of  absorption  of  sound 
in  the  water,  attenuation  by  the  rough  ice  boundaries,  and  attenuation  of 
sound  by  the  bottom. 

Two  computer  programs  were  written  in  Fortran  IV  to  evaluate  the  liquid- 
and  solid-bottom  integral  solutions  derived  by  Kutschale  (1970,1972).  A 
comparison  of  the  FFP  computations  with  those  computations  by  normal  mode 
theory  from  the  corresponding  integral  solutions  are  nearly  identical,  but 
the  FFP  technique  is  far  more  convenient  and  at  least  an  order  of  magnitude 
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faster  since  the  computations  are  done  directly  from  the  integral  solution 
without  first  computing  the  roots  of  the  frequency  (period)  equation  and 
then  summing  the  normal  modes. 

A  number  of  comparisons  of  propagation  loss  computations  have  been  made 
using  ray  theory  (the  program  developed  by  Aerophysics  Research  Corpora¬ 
tion)  ,  normal  mode  theory,  and  the  numerical  integration  method  developed 
by  Lamont- Doherty.  On  page  89  are  shown  the  computed  propagation  loss 

curves  against  range  by  various  methods.  Also  included  are  the  loss 
curves  of  Mellen  and  Marsh  and  of  Buck  (see  References  4  and  12  of  Section 
4).  These  results  indicate  that  for  long  range  propagations  the  losses 
calculated  by  ray  theory  agree  well  with  those  by  other  methods. 

Details  of  the  numerical  integration  method  are  summarized  in  a  separate 
report  (Reference  5). 
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6.  COMPUTED  RESULTS  AND  PLOTS 


In  this  section  we  present  the  final  computed  results  in  graph  form.  The 
detailed  output  is  reproduced  in  Part  II  of  the  report. 

As  has  been  mentioned  elsewhere  in  the  report,  the  frequency  of  propagation 
is  20  hz.  The  location  of  sources  and  receivers  is  giver,  and  plotted  on 
page  57  -  page  76.  The  typical  temperature  and  salinity  data  are  the 
recent  measurement  values  taken  at  Base  Camp  Station  (72°-02.0N  latitude 
and  211°-24.8  longitude): 


Depth  (meters) 
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TABLE  OF  SOURCE  AND  RECEIVER  LOCATIONS  AND  BOTTOM  DEPTHS 
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CALCULATED  RESULTS  AND  BEST  FIT  CURVE  FROM  .MEASUREMENTS  (B.  M.  BUCK) 
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